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Abstract 

A challenge for nuclear physics is to measure masses of exotic nuclei up to the 
limits of nuclear existence which are characterized by low production cross 
sections and short half-lives. The large acceptance Collector Ring (CR) [l[ at 
FAIR 0] tuned in the isochronous ion-optical mode offers unique possibilities 
for measuring short-lived and very exotic nuclides. However, in a ring de- 
signed for maximal acceptance, many factors limit the resolution. One point 
is a limit in time resolution inversely proportional to the transverse emit- 
tance. But most of the time aberrations can be corrected and others become 
small for large number of turns. We show the relations of the time correction 
to the corresponding transverse focusing and that the main correction for 
large emittance corresponds directly to the chromaticity correction for trans- 
verse focusing of the beam. With the help of Monte-Carlo simulations for the 
full acceptance we demonstrate how to correct the revolution times so that in 
principle resolutions of Am/m = 10^^ can be achieved. In these calculations 
the influence of magnet inhomogeneities and extended fringe fields are con- 
sidered and a calibration scheme also for ions with different mass-to-charge 
ratio is presented. 
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1. Introduction 

The Collector Ring of the FAIR facility is a symmetric, achromatic ring 
with two arcs, two straight sections and a total circumference of 221.5 meters. 
It is designed for operation at a maximum magnetic rigidity of Bp = 13 Tm 
[l|. It will be operated in three ion-optical modes, two of them providing 



Preprint submitted to Nuclear Instruments and Methods A 



May 27, 2013 



fast pre-cooling of either antiprotons or radioactive ion beams In the 
third mode (isochronous optics) the CR will be operated as a Time-of-Flight 
(ToF) spectrometer for the mass measurement of exotic very short-lived nu- 
clei (Ti/2 > 20/is) produced and selected in flight with the Super-FRS frag- 
ment separator [Sj. This technique for mass measurements has been devel- 
oped at the ESR at GSI j3] and meanwhile was also used successfully in the 
CSRe at IMS Lanzhou and the RIKEN Rare-RI ring is under construction 

The first basic condition for isochronicity is the well known equation j3] 
which describes the revolution time (T) of ions stored in the ring and its 
deviation (AT) depending on the difference in mass-to-charge ratio {m/q) or 
different velocity (v): 



AT 1 A{m/q) fY aAv AT, 



+ (^-1)— + ^^, (1) 



T 7^ (m/q) \jf J v T 

where 7 is the relativistic Lorentz factor and 7^ is the transition energy of 
the ring. The isochronous condition is reached when 7 = 7^. This means, 
the second term in Eq. ([T]) vanishes and A(m/g) can be determined from 
the observed AT. It is remarkable that T is directly proportional to m/q. 
But in addition to the intended time differences for different masses and the 
important first order isochronicity condition the resolution depends on more 
additional influences on time-of- flight which are collected in ATothers- These 
aberrations can also have higher order contributions. Effects of nonlinear 
field errors, fringe fields of magnets, closed orbit distortion and the transverse 
emittance negatively act on the resolution [7]. 

The general layout of the CR has been described in ref . Q . The isochronous 
mode has been calculated for different values of 7f pl. In this paper we will 
always refer to the mode for 7^ = 1.67 which allows to measure masses 
up to m/q = 3.1u/e. In this case the momentum acceptance becomes 
Ap/p = ±0.5% when having at the same time a transverse acceptance of 
100 mm mrad in both planes. 

A correction of time-of-flight aberrations due to a large transverse phase 
space distribution is essential or otherwise the goal of large phase acceptance 
in the CR cannot be achieved with sufficient time resolution [3]. The simple 
reason for this influence is that ions which oscillate around the optical axis 
have a longer path length for one turn which leads to a corresponding time 
spread. As it was indicated before [lo| and will be also shown in this article 
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these time-of-flight aberrations are directly related to the chromatic aberra- 
tions of the optical system. It is known that in a straight system with only 
magnetic lenses defined by magnets outside of the beamline a correction is 
not possible [ll|. With additional electrostatic elements this would be pos- 



sible [12|, however, these are not available for the high rigidity of the CR 



beam. At least a correction in a system with a bent optical axis is possible 



12[. A recipe for such a full second or higher order achromat was given by 



Brown [13||, and partially this was also applied in the spectrometer TOFI 



14j | for time-of-fiight mass measurements of nuclear masses. However, when 



a storage ring is constructed in this way a four fold symmetry system which 
is imaging violates the basic stability criterion of a nonzero phase advance 



for a closed ring [15 



This is also the result of a general approach to higher order achromats by 



Wan [16[, the minimum symmetry needed for a second order achromat is a 



four-fold symmetry, but based on the Lie algebra description of the system 



17l | it is found that the tunes of the individual cells should be integer. A good 



approximation to a second order achromat may be found with the method 
of extending a full achromat by straight matching sections for adjusting the 
tunes. In this case the overall second order geometric aberrations can be kept 
zero as pure quadrupoles are free of them and by an overcorrection with the 
sextupoles at least the tune shift be corrected (so-called pseudo achromat) 



18[. As the matching section can be short compared to the whole arc section 
this usually leads to systems close to a full achromat. 

This dilemma may not be so crucial in a ring for only few revolutions 
but for commissioning and other detection techniques like Schottky pickups 
stability of motion must be guaranteed. In this paper we will show a solution 
that does not lead to full achromaticity over one turn but to a correction in 
the limit of many turns. By averaging the time of fiight over many turns still 
a good resolution for a large phase space can be obtained. This less stringent 
requirement also helps to avoid the usual problem of large higher order field 
components, which in turn could impose other higher order aberrations. 

In the CR we want to use ToF detectors inside the ring which do measure 
the revolution time for each turn. This means we must observe the ions over 
a series of turns and see how large the remaining error will be. This is anyway 
necessary as the resolution of the detector itself is limited and only averaging 
over a larger number of turns helps to obtain the desired resolution. 
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2. Calculation of Second-Order Time Deviations 



The relative revolution time difference between an arbitrary and the ref- 
erence particle can be expressed in terms of the initial coordinates as a Taylor 



series in a second-order approximation |19l . |20 
AT T - Tf 



{t\x)cX + {t\a)ca + {t\6)cS + {t\xx)cX 



T To 

+ {t\xa)cxa + {t\aa)ca^ + {t\yy)cy'^ + {t\yb)cyb 

+ {t\bb)cb^ + {t\x6)cx6 + {t\a6)ca6 + {t\SS)J^, (2) 

where {x),{y) are the transverse position coordinates and (a), (6) the trans- 
verse momenta divided by the forward momentum of a reference particle (po)- 
The coefficients of the Taylor series correspond to the partial derivatives of 
the first coordinate before the line with respect to the following. The index c 
marks coefficients normalized by the total time-of-fiight nTo, where n is the 
number of turns and Tq is the revolution time of the reference particle. The 
fractional momentum deviation 6 is given hj p = Pq{1 + 6). 

In the first-order achromatic ring, with {x\6) = {a\S) = 0, the first-order 
transverse time coefficients {t\x) and {t\a) vanish simultaneously [lo[. The 
necessary condition to be an isochronous ring in first order is {t\S) = 0, i.e. 
7 = 7j. The second-order isochronous condition is fulfilled when {t\66) = 0, 
which can be corrected with one family of sextupole magnets installed in a 
dispersive section of the ring. 

With such a good adjustment of 7t even to higher order in an achromatic 
ring the largest contribution to AT comes from the second-order geometric 
time aberrations. For a beam of one species in an ideal isochronous ring 
without higher order field errors or closed orbit distortions only pure betatron 
motion exists. In such a ring the time spread is directly related to the 
transverse emittances (^x,?/) [2|: 

(^)_^^ -^(..<7.>+..<7.», (3) 

where ex,y is the transverse beam emittance, which can be regarded as the 
ring acceptance, and < 7a;,y > are the Twiss parameters '-yx,y averaged over 
the whole circumference of the ring. This result simply describes the longer 
path length of ions traveling with larger betatron amplitude over one turn. 

From the maximum deviation given in Eqs. [H [3] one can derive the mass 
resolution over the whole beam emittance. For the CR, with acceptance 
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100 mm mrad in both planes, the hmit of mass resolution would be about 
10~^, which is insufficient for precise mass measurements. Therefore, in order 
to reach the necessary resolution of 10~^ the transverse emittance would have 
to be limited to 10 mm mrad in both planes. As a result, the transmission 
of the ions into the ring would be reduced drastically. However, the mass 
resolving power can be improved using second-order corrections and keeping 
the transverse emittance large. 

For the second order terms it is useful to investigate their behavior over 
many revolutions. As the system is bound and the betatron oscillations 
have a non-zero phase advance the initial coordinates of a single ion change 
from turn to turn within fixed limits. The time deviations over many turns 
simply sum up. The fluctuations are periodic with the tune and have an 
average value of zero. Therefore, the time deviations given by a single matrix 
coefficients which depend on mixed initial coordinates will oscillate with a 
fixed maximum amplitude. When looking at the relative time difference 
over many turns these contributions will become small. The quadratic terms 
however will continue to increase from turn to turn and the relative time 
deviation converges towards a constant value. FigJT] illustrates this behavior. 
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Figure 1: Evolution of the relative second-order geometric time aberrations as a function 
of the number of turns in the CR for a single ion with start positions x and a on the edge 
of the acceptance of 100 mm mrad. The inserted numbers show the limit over many turns. 
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The constant second-order geometric aberrations in time of EqJ2]are most 
significant. Without correction their contribution corresponds to the right 
part of Eqj3] proportional to the transverse emittance. 



3. Correction of Geometric Time Aberrations and Chromaticity 

Let us assume a beam of one species circulating in the ring turn by turn. 
We observe it in the symmetry plane of the ring where the phase-space ellipse 
is upright (Twiss parameter a = 0) and this condition is restored after each 
turn. For simplicity we inspect only the horizontal plane, the arguments for 
the vertical plane are the same. The transfer matrix for the full circumference 
can be expressed by the betatron function /3 and the phase advance /i: 

/ cos/i (3 sin fi\ _ f {x\x) (x|a)\ , . 



y -gSm/i cosyu J \\Ci\x) \a\a)^ 

The matrix elements {t\xx), {t\xa) and {t\aa) can be expressed for an 
achromatic ring as 20|, |21| : 

(5) 
(6) 
(7) 

where A is a constant with the dimension of inverse velocity. The coefficients 
in Eqs. (jSHZ]) correspond to the terms of the transfer matrix M. 
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The terms {x\a5) and {a\x5) are connected via 21 



{x\a5) (a|x) = {x\x) {a\a5) — {a\x5) {x\a) + {x\x5) (a|a). (8) 
Combining Eqs. El El El we obtain: 

it\aa) +/3^(t|xx) = -^[(x|x(5) + {a\a6)]. (9) 

The shift of the tune (Qox) of the ring in relative units for a certain 
momentum shift 5 is called relative chromaticity {^ix)- It can be written as 
(see appendix): 

= . ■ [{x\x5) + {a\a5)]. (10) 
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Therefore, by combining Eqs. ([21 ITU]) one obtains: 



{t\aa) + (3^ {t\xx) = -47rA/3Qox6: 



(11) 



The emittance of the beam is 



(12) 



and then Eq|TT] can be written as: 



(t\aa) + {t\xx) xl 



(13) 



As one can see the remaining time error for an ion starting in both position 
and angle exactly with the coordinates like the half axis of a matched ellipse 
{xq and Oo) is still proportional to the emittance but also proportional to the 
chromaticity which can be corrected to zero (^i^. = 0). 

Next we want to show that for all other starting points in phase space at 
least the relative time deviation over many turns becomes zero. In general, in 
a first order achromatic ring with mirror symmetry the two chromatic matrix 
coefficients coincide, {x\x5) = {a\a5) |22|. This means the chromaticity 
correction directly demands, see EqJTOl 



The time deviation by the mixed term {t\xa)c converges to zero for large 
number of turns. 



When inserting this into the symplectic relation (Eqj6]), one sees that also 




(14) 




(15) 








(16) 



and by inserting this into EqJT] also 





From this immediately follows also 




(18) 
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This means the correction of chromaticity corresponds directly to the 
correction of the otherwise non- vanishing time deviations. For this task only 
one condition has to be fulfilled and one family of sextupole magnets is 
needed. Together with the same condition for the vertical plane, {y\y5) = 
0, and the necessary (t\66) = 0, this requires in total three independent 
sextupole correctors. Then we can reach a regime where all time deviations 
vanish in the limit of many turns and the isochronous ring behaves like a 
second-order achromatic system in relative time deviation. 



4. Further Correction of Dispersion 

While the first-order achromaticity of the ring is a basic prerequisite, there 
are also higher order contributions to dispersion. The mixed aberrations 
{t\x6)cx6 and {t\a6)caS do vanish in the limit of many turns, but still for 
a short observation time one wants to keep their magnitude small. Again 
we establish a relationship with the corresponding chromatic aberrations in 



transverse space. The symplectic relations for the achromatic ring yield [20 



21 



{t\x5) = X ■ [{x\x){a\55) - {a\x){x\55)] , (19) 
{t\a6) = X ■ [{x\a){a\66) - {a\a){x\66)] , (20) 



and 22 



ia\SS){x\a) - {x\5S) [1 - {x\x)] = (21) 

In a ring with non-zero phase advance the first order coefficients in Eqj2T] 
cannot be zero, and the only solution is 

{x\66) = {a\66) = , (22) 

which then over a single turn leads directly to 

{t\xS) = , (23) 
{t\a6) = . (24) 

This can be achieved with one additional sextupole which corrects the 
nonlinearity of dispersion after one turn. The influence will be demonstrated 
in section 6 by a Monte-Carlo simulation. Fig|2] gives an overview on one 
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Figure 2: Layout of one quarter of the CR. The different magnet types are indicated and 
the numbers mark the sextupole magnets which are used in the correction schemes with 
only one, three or four sextupoles. Inside the apertures dispersion curves are shown. In the 
fuU ring with four-fold symmetry each multipole magnet belongs to a family of magnets 
with same value. 

quarter of the CR. The positions of sextupoles used in different correction 
schemes are indicated, as well as additional superimposed octupole coils in- 
side some of the quadrupole magnets. 

We have seen that after one turn the ring must be achromatic at best 
also to higher orders. But on the other hand in the arcs we need a well 
defined dispersion. This must even be a bit non-linear because the Lorentz 
factor 7 is not constant. It does not even depend linearly on S, and we want 
the condition 7 = 7t (see Eq{T]) to be fulfilled over the whole momentum 
acceptance. The nonlinear dispersion then defines the small but necessary 
term {t\66). In practice this demands symmetric sextupoles for correction to 
first create and to adjust dispersion in the arcs and then to compensate it 
again. 

In the non- dispersive vertical direction the chromatic terms described 
above do not exist in a second-order approximation. What remains are the 
purely geometric terms not in diagonal direction inside the phase space el- 
lipse, {t\xa) and {t\yb), for which also a relation to the chromatic aberrations 
in positions exists, see EqEl But as we cannot have directly a full second- 
order achromat, they are only corrected in the limit of many revolutions. 
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5. Field Quality 



5.1. Dipole Homogeneity 

The dipole field homogeneity is crucial for the time resolution. It can be 
described by a series of higher order field indices rii 



The sector radius po is 8.125 m and the given values refer to the case of 
maximum field strength {B = 1.6 T). The chromaticity depending only 
on the quadrupole fields in the ring, usually called "natural chromaticity", 
is one contribution, but the magnet design predicts a large second order 
deviation, which causes an effect even stronger than the natural chromaticity, 
making chromaticity worse in horizontal direction and overcompensating it 
in vertical direction. But as n2 represents a pure sextupole component it can 
be completely corrected with three sextupoles including the {t\66) term and 
with no increase in maximal sextupole magnet strength compared to the case 
of a homogeneous dipole field. 

By exciting one assumed multipole component of order i in the main field 
of all dipoles one sees that each multipole component leads to a corresponding 
time error {t\S^) after one turn. These deviations add up over many turns and 
do not average out. Compensation must be made with adjustable multipoles. 
So far sextupoles and octupoles superimposed to some of the quadrupoles and 
one decapole are foreseen, see Fig|2J 

Figl3] shows the remaining field shape deviations of order four and higher. 
For demonstration the n2 component in the dipole was assumed to be com- 
pensated so that only higher order terms (mainly n^j contribute. Sextupoles 
were adjusted to overcompensate the second-order dependence to minimize 
the overall field deviation. Still the time spread assuming a momentum distri- 
bution according to the full CR acceptance would be limited to (Tt/T ^ 10~^. 
Clearly, a weak decapole field close to the maximum of dispersion makes com- 
pensation more effective and easier to adjust. 

In a simplified model the relation between the relative field quality AB / B{x) 
and the time deviations AT/T{6) can be described with an effective disper- 
sion of (x|5)eff = 29 m. Then the ratio becomes 



B{x) 




(25) 
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Figure 3: Comparison of dipole inhomogeneity from Eg 1251 as function of position x with 
the resulting time deviation as function of 5. For comparison the curves are shifted on top 
of each other with the help of an effective dispersion coefficient (a;|(5)cff = 29 m. 



This means the required field homogeneity must be very close to the desired 
mass resolution. As a dipole magnet with large gap will not provide this 
many multipole correctors must be foreseen. 



5.2. Fringe Fields 

The fringe fields of the magnets include also many more multipole terms 
and the CR is designed for large apertures. Therefore, the infiuence of the 
fringe fields was studied in detail. The analysis is based on a realistic fringe 
field distribution as calculated for the magnet design of the CR dipole and 



quadrupole magnets [23(]. They were approximated by Enge functions and 



corresponding transfer matrices were calculated using the program of ref. 24 



An alternative is the method of fringe field integrals 25|, l26|] which yields 
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ractically the same results. Both are implemented into the GICOSY code 



[2 71 . |28| used for the analysis. 

The simple quasi- analytic fringe field integral equations reveal that the 
first influence of the quadrupole fringe fields lies in a change of position and 
angle when describing the effect by a homogeneous main field with an addi- 
tional short fringe field matrix. Each entrance and exit matrix contains terms 
{x\x) = — (a|a) = —{y\y) = {b\b) = 1 ± w with w being proportional to the 
gradient and aperture radius of the magnet. In case of the large quadrupoles 
of the CR with aperture radius tq = 0.168 m and the strongest excitation 
{B = 0.4898 T for Bp = 13 Tm) the value can reach w = 0.000727. Over 
the whole ring this still leads to a matched beam with slightly altered beam 
parameters and only small dispersion mismatch. When including them the 
tunes of Qx = 2.153 and Qy = 3.186 shift by —3.7% or —1.2%, respectively. 
Even without a readjustment of the quadrupole setting the compensation for 
different velocities {t\6) still stays on the level of few times 10~^. Second order 
terms are not significantly affected, only slightly by the natural chromaticity 
of the changed first order focusing, see tabl^ 

However, the well known third order aperture aberrations of the quadrupole 
fringe fields make a significant difference. After the described compensation 
of the second-order path length differences and counterbalancing of dipole 
inhomogeneities with other multipoles, which is largely possible for the given 
numbers, they introduce the strongest limitation of ToF resolution. The cal- 
culations show that this is largely independent of the gap size and fringe field 
shape. The main contributions exist even for small aperture radius with a 
very short extension of the fringe field as long as the shape follows a realistic 
shape and not a simple cut off which actually violates Maxwell equations. 

The dipole magnets in CR have no edge angle for additional focusing. 
By the nature of the fringe fields some additional focusing in first order in 



vertical direction is added [25[. But its magnitude is smaller than in the 
quadrupole fringe fields. However, in any magnetic sector field the fringe 
field introduces an additional focusing in vertical direction dependent on the 
horizontal beam direction. Its size is practically constant as in leading order 
it directly depends only on the deflection radius, {b\ya) = l/po (isf. This 
couples the vertical focusing to the horizontal dispersion and changes the 
vertical chromaticity. 

In addition the vertical focusing itself is also associated with chromatic 
terms depending on the shape of the field distribution. When adding realistic 
dipole fringe fields to a CR with adjusted first order tunes they shift the tune 
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Table 1: Chromatic tune shifts for S ~ 0.5% by different influences separately. The shifts 
are given relative to a ring with no chromaticity. 







AQy 


natural chromaticity (sext. off) 


-0.0196 


-0.0138 


dipole field with n2=1.15 


-0.0260 


0.0176 


dipole fringe fields 


< 10-4 


0.0032 


quadrupole fringe fields 


< 10-4 


< 10-4 



in vertical direction by up to 

^Qy = Qoy ^max ■ (28) 

This shift is small compared to the natural chromaticity, see table [1], and 
the tune in horizontal direction is practically not affected. The sign of the 
chromaticity from the dipole fringe fields even compensates partially the 
natural chromaticity. All which is required is a smaller readjustment of the 
sextupole setting. 



In general most sensitive to higher order multipole components are the 
terms for isochronicity depending purely on 6. This means after each ad- 
justment of chromaticity to reduce the transverse path length differences one 
sextupole, octupole and decapole must be adjusted. The size of the other re- 
maining chromatic transverse path length differences in third or higher order 
can amount to AT /T ^ 10-^ depending on the exact multipole correction. 

6. Results of Numerical Simulations 

Numerical simulations were done by tracing single ions directly inside 
the GICOSY program by transforming their coordinates by up to fifth order 
transfer matrices for each optical element. This method differs from inserting 
the initial coordinates in one overall transfer matrix by not truncating higher 
orders which appear in a multiplication of lower order matrices with order 
larger than one. In the simulations an initial five- dimensional phase space 
area was assumed with an emittance of 100 mm mrad in each plane and a 
momentum spread of ±0.5%. The shape of the initial ellipses was chosen to 
have matched beams in x and y plane. The ellipses in both planes were filled 
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independently. The stepwise tracing of ions also makes insertion of collima- 
tors possible, however for the chosen size of phase space the transmission is 
larger than 98% and in the results shown no collimators were used. This 
approach is similar to the program MO CADI [29,] which was used as a cross 
check up to third order. 

The averaging of time aberrations over the full phase space allows a real- 
istic comparison of the time resolution. To better understand the resolution 
over a lower number of turns Figj4] shows the standard deviation in time 
as function of turn number in a simplified ring with perfectly homogeneous 
dipole magnets and sharp cut off multipole fields. As one can see, conver- 
gence to an almost constant limit goes fast, which justifies our approach of 
correcting mainly terms not vanishing over many turns. When correcting 
only the velocity dependence the time resolution for large emittance remains 
still too large for good mass measurements of heavy ions. The chromatic- 
ity correction improves this limit by almost a factor of ten. With three 
sextupoles there are small oscillations in the time spread which result from 
non-matched dispersion in higher order with the corresponding time errors, 
see EqJ2^ The size of this dispersion oscillates with the frequency of the 
horizontal tune. A fourth sextupole can correct this and thereby make con- 
vergence faster. However, more and stronger sextupoles are needed which 
inflict more higher order aberrations. What remains after many turns are 
time errors mainly of higher order, even though the simulations already in- 
clude one octupole family located near the maximum of dispersion in the arc 
of the ring. It is needed to correct the largest third-order contributions to 
the revolution time, which is the purely chromatic term {t\666). 

The expected imperfections of the dipole magnets require stronger sex- 
tupole corrections than without and also a readjustment of at least one oc- 
tupole for the condition {t\666) = 0. Then the decapole component of the 
dipole dominates. Already Fig J3] makes it clear that without adding one sepa- 
rate weak decapole the goal in resolution cannot be reached. The remaining 
ToF deviations for different field imperfections are shown in Figj5] by the 
standard deviation in ToF as function of momentum deviation. In these 
calculations a correction of chromaticity with sextupoles and a correction of 
the pure momentum dependence up to fourth order was considered with the 
help of three sextupoles, one octupole and one decapole. As fringe fields with 
real extension were inserted also a small readjustment of quadrupole and sex- 
tupoles was required. The influence of the additional multipoles in the dipole 
magnet can mostly be corrected but in combination with the large emittance 
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Figure 4: Standard deviation of relative ToF as function of turns. The Monte-Carlo 
simulations show three cases: one sextupole to correct only {t\66), three sextupoles to 
correct in addition chromaticity in x and y plane, and 4 sextupoles to correct as well second 
order dispersion {x\SS). In one calculation (green curve) fringe fields on the quadrupoles 
were considered in addition. 

some contributions remain which depend on the actual sextupole setting. 
Here values for the simple case with only three sextupoles are shown. Also 
the fringe fields limit the ToF resolution. The influence of the quadrupole 
fringe fields dominates. In general limitations from higher order mixed com- 
ponents remain. The dependence on transverse motion alone is corrected 
very well by sextupoles but in combination with dispersion more deviations 
show up. 

That quadrupole fringe fields have a basic unavoidable influence is demon- 
strated by adding them to the perfect pure multipole fields in FigJH To 
show the basic influence the pole tip radius size was reduced very much so 
that the fringe field extension becomes very short while keeping the gradient 
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Figure 5: Relative time deviation after 100 turns for the full accepted transverse phase 
space as function of momentum deviation S. After correction of pure S dependence and 
chromaticity with sextupoles the remaining influence of including dipole inhomogeneities 
and quadrupole or dipole fringe fields is shown. The green curve shows all effects combined. 

constant. In this case the first-order focusing becomes identical to the case 
without fringe fields. However, in third order many aberrations remain which 
do not vanish for a real field distribution and are even present for very small 
gap size. As apertures were not used in the calculation the CR with small 
quadrupole radius still has the same acceptance. One can see that even then 
the quadrupole fringe field represents a major contribution to the limit in 
ToF resolution. 

Considering all influences the dependence on turn number is shown again 
in Figjni now with all fringe fields and dipole inhomogeneities. Applying in 
addition further octupole corrections also some of the third order ToF de- 
viations can be corrected but the effect of quadrupole fringe fields is more 
difficult to compensate. Contrary to the dipole inhomogeneity it does not 
correspond to a simple multipole component which can be compensated with 
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an opposite field nearby. More terms coupling the longitudinal and trans- 
verse motion become important. With additional four weak superimposed 
octupoles the result is improved a bit, see Fig. [61 One can reach a resolution 
of up to ar/T ^ 3.0 ■ 10~^, which corresponds to a mass resolving power of 
iTT'/o'm = 1-2 ■ 10^ or o"m = 100 keV at m = 130 m. Further improvement in 
mass precision can be reached by detecting many ions of one species. 




0.0 H 1 1 1 1 

5 10 15 20 

turn 



Figure 6: Time spread as function of turn number, similar to Fig|4]but with realistic dipole 
inhomogeneities and fringe fields. In addition to only one octupole also four octupoles 
families were used for correction. 



7. Calibration and Tuning 

In principle ions with different m/q can be separated if their mean time 
separation AT is larger than the full time spread of the beam over a number 
of turns and including detector resolution. However, only ions of very similar 
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m/q can be in a good velocity region to be isochronous. For other m/q but 
same magnetic rigidity Bp the velocity 7 must be different from 7^. In prac- 
tice this leads to the problem of finding enough calibrants with well known 
masses. Assuming perfect isochronicity for an ion of m/q = 140/50 u/e leads 
for m/q = 139/50 u/e to a large time spread of ar/T = 1.1 ■ 10~^ in case of 
full acceptance of the CR with as = 0.003. 

At the same time a method for tuning of isochronicity with all the men- 
tioned corrections must be found. Scanning over the momentum acceptance 
with a narrow emittance beam of variable energy is used in the ESR but 
there is no electron cooler foreseen in the CR and stochastic cooling will not 
work in the isochronous mode. Instead a second ToF detector on the straight 



section of the CR shall be installed 30 . From the time difference between 



the two detectors the velocity can be deduced, whereas over a full turn the 
revolution time does not depend on velocity. From the velocity the magnetic 
rigidity can be calculated and by following many ions of different velocity a 
full isochronicity curve can be recorded. In addition with the help of slits 
and steerers in the preceding Super-FRS separator ^ and transfer beamline 
also regions of the transverse phase space can be selected. 

The required resolution in time-of-fiight for this velocity measurement is 
not as high as for the exact mass determination. It is sufficient to provide 
enough calibration points in 6 over the whole CR momentum acceptance 
and to be able to distinguish ions of different mass number. Even with 
masses only from theoretical predictions the uncertainty in the conversion 
from velocity to Bp and 6 is small enough. 

At best the detectors are located on two ends of a straight section. The 
distance in case of the CR is 22 m. The time deviation of the present ToF 



detector was determined experimentally and amounts to at = 37 ps [31[. This 
leads to a time resolution of at/t = 6 ■ 10~^, which is improved by observing 
an ion over many turns. Even when considering the limited efficiency of the 
detector, which will not give a valid time signal for each turn, this will still 
yield At/t > 3 ■ 10~^ after 50 revolutions. At 7 = 1.67 this corresponds to 
an error in 5 of 7 ■ 10^^. Compared to the ring acceptance of ±5 ■ 10~^ this 
gives a good resolution. On this level the described second order path length 
differences which still exist over the straight section are not critical. 

It is noteworthy that the energy loss of ^^''Sn ions in two foils of 20 pg/cm^ 
over 50 turns corresponds to 5 = —8.6 ■ 10^^. That means it is possible to 
trace the energy loss in the foils. One could even evaluate the revolution time 
only in a narrowed region that corresponds to good isochronicity. By this 
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Bp tagging also ions with larger difference in m/g could be included in the 
analysis. One would simply evaluate them only when they are close to 5 = 0. 
By the energy loss in the foil even ions with initially too high Bp will pass 
this point if the difference is less than the CR acceptance. This would lead 
to a large gain in acceptance compared to Bp tagging as performed today by 



selection with slits inside the FRS or ESR 32 



8. Conclusion 

We have presented a method for correcting second order time-of-flight 
aberrations in a storage ring. As a full correction over one turn is not feasible, 
instead a correction over many turns is proposed. The necessary requirements 
on the higher-order multipole correctors were derived. A Monte-Carlo simu- 
lation for the full phase-space acceptance shows that the main part of the not 
corrected terms average out fast. These correction must take into account 
the unavoidable influence of dipole inhomogeneities and of extended fringing 
fields which were investigated in detail for the CR. Other multipole com- 
ponents can be compensated with more multipole correctors but especially 
the spherical aberrations in the fringe field of quadrupoles are hard to com- 
pensate. Nevertheless, the achievable time resolution is sufficient for mass 
measurements even for the full acceptance of the CR. A detection scheme 
with two detectors in the ring helps in tuning and by tagging properties of 
individual ions the accuracy can be improved. 



9. Appendix 

Here we derive the expression for the nonlinear chromaticity described in 
matrix coefficients. The phase advance can be expanded with respect to 5 
as 
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= XI cos = X Xn5"- (29) 

n=0 n=0 

Since cos p, is the trace of the transfer matrix M (see Eq. H]), which can also 
be expanded with respect to 5 one obtains 

X„ = ^Tr(M„), (30) 
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where M„ is in terms of the transfer matrix coefficients: 

M»=f HI"! <!K:>Vn = 0. 1. 2. 3... (31) 



(a|x(5") (a|a(5") 

Thus, from Eqs. fl29|3ip one can extract the phase advance described in matrix 
coefficients and correspondingly obtain the chromaticity: 
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